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ABSTRACT 

We analyze photometric data in SDSS-DR7 to infer statistical properties of faint satellites associated 
to isolated bright galaxies (M r < —20.5) in the redshift range 0.03 < z < 0.1. The mean projected 
radial number density profile shows an excess of companions in the photometric sample around the 
primaries, with approximately a power law shape that extends up to ~ 700 kpc. Given this overdensity 
signal, a suitable background subtraction method is used to study the statistical properties of the 
population of bound satellites, down to magnitude M r = —14.5, in the projected radial distance range 
100 < r p /kpc < 3 < R v i r >. The maximum projected distance corresponds is in the range 470 — 660 
kpc for the different samples. We have also considered a colour cut consistent with the observed colours 
of spectroscopic satellites in nearby galaxies so that distant redshifted galaxies do not dominate the 
statistics. We have tested the implementation of this background subtraction procedure using a 
mock catalogue derived from the Millenium simulation SAM galaxy catalogue based on a ACDM 
model. We find that the method is effective in reproducing the true projected radial satellite number 
density profile and luminosity distributions, providing confidence in the results derived from SDSS 
data. We find that the spatial extent of satellite systems is larger for bright, red primaries. Also, 
we find a larger spatial distribution of blue satellites. For the different samples analyzed, we derive 
the average number of satellites and their luminosity distributions down to M r = —14.5. The mean 
number of satellites depends very strongly on host luminosity. Bright primaries (M r < —21.5) host on 
average ~ 6 satellites with M r < —14.5. This number is reduced for primaries with lower luminosities 
(—21.5 < M r < —20.5) which have less than 1 satellites per host. We provide Schechter function 
fits to the luminosity distributions of satellite galaxies where the resulting faint end slopes equal to 
1.3 ± 0.2, consistent with the universal value. This shows that satellites of bright primaries lack an 
excess population of faint objects, in agreement with the results in the Milky Way and nearby galaxies. 



1. INTRODUCTION 

According to the currently accepted model of struc- 
ture formation, gala xy systems arise as the result of hi- 
erarchical clustering (|White &: Fr enk 1991]; iBertsc hingcr 
11994 iCole et al.l [199J). The details by which galaxies 
form and evolve in dense or moderately dense environ- 
ments, where galaxy- galaxy interactions are frequent and 
matter distributes in a rich substructure, depend on the 
characteristics of those environments. The assembly of 
galaxy systems entail the process of matter accretion, 
governed by gravity, as well as astrophysical phenom- 
ena, such as the efficiency of gas to cool and collapse 
or the energy feedback related to the late stages of stel- 
lar evolution (jViola et al.ll2008HKang et al.ll2006l ). Some 
of these details are not yet fully understood, and obser- 
vational evidence is fundamental to constraint structure 
formation and evolution models, specially because faint 
galaxies are more sensitive to astrophysical processes like 
supernova feedback and ram pressure stripping. In par- 
ticular, statistical studies of systems of galaxies are key 
to understand the transformations of galaxies due to 
the interactions be tween galaxies and their environment 
(|Nichol et all 120031 ) . 

Although the formation of a galaxy is believed to take 
place on the potential well of a dark matte r halo, not all 
haloes host a galaxy. As pointed out first by lKlvpin et al.l 
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(1999), the number of halos in numerical simulations 
were an order of magnitude greater than the observed 
number of satellites in the local group dMoore et al.lll999l : 
iKravtsovetafl l200l iStrigari et al.l 120071 ). This differ- 
ence has been a matter of lively debate, being attributed 
either to a lack of observed g alaxies or to an exce ss of 
formed objects in simulations. iWillman et al.l (|2002f) dis- 
cussed the possibility of under counting satellite galaxies 
in the Milky Way and estimated the actual number of 
satellites in about twice the known population at that 
time. The authors argued that galactic extinction and 
stellar foreground can lead up to 33 per cent of incom- 
pleteness, and then the number of Milky Way satellites 
at low galactic latitude an d at galacto-ce n tric d istance 
might be underestimated. iSimon fc Gehal (|2007l ) found 
that the number of satellites in the Milky Way system 
was greater than previously expected, based on an anal- 
ysis of the SDSS data. With these findings, the dis- 
crepancy between the number of observed and simulated 
satellites reduces to a factor of nearly 4. Nevertheless, 
applying a background subtr action technique on photo- 
metric data from SDSS-DR7, iLiu et al.l (|2010t ) found re- 
cently that the Milky Way galaxy has significantly more 
bright satellites than a typical galaxy of its luminosity. 

The study of the spatial distribution of satellites 
around primaries and clusters has become favored by in- 
creasingly large galaxy redshift surveys. Many works ad- 
dress observational studies of radial distribution of galax- 
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ies in spectroscopic samples dCoil et al J 120061: iLin et al.l 
[2001 lYang et alJl2005t ICollister fc Lahavl 1200511 . mostly 
around galaxy clusters , and on bright p rimary galaxies 
([Sales fe L amb as 2 005t iChen et alJ l2006) . Deeper sam- 
ples have been also used to compute pro jected density 
profil es, using background subtraction ([Hansen et al.1 
2005!) around MaxBCG galaxy systems. Moreover, 
galaxy projected density profiles were computed based 
on projected correlatio n function det erminations in red- 
shift galaxy cata logues dLi et al.f 2007) and deeper photo- 
metric samples ([Wang et alJl2010l ) around a set of spec- 
troscopically identified galaxies. 

The distribution of satellite luminosities is also key to 
the development of models a nd understanding of the pro- 
cesse s of galaxy formation ([Benson et al.l 120031 : iBensonl 
120101 lOkamoto etall 12010( 1 . Given the low luminos- 
ity of most satellites, however, their observation is usu- 
ally one rous an d hard ly accessible outward of the Local 
Group. iMateol ([1998H put forward a detailed census of 
dwarf galaxies, from which a flat faint end of the luminos- 
ity function in the Local Group could not be discarded. 
Background subtraction methods have been widely used 
to obtain galaxy luminosity function in clusters, and 
also results of this procedure on individual clusters have 
been reported (e.g. lOemlerlll974D . Since it is not limited 
to the computation of luminosities, it can be also use d 
to obtain the c o lor-m agnitude relation ([Pimbbleti [20081. 
lAndreon et al.l ([2005(1 present a variation of the back- 
ground decontamination method, avoiding the use of ar- 
bitrary binning and incorporating the background noise 
as part of a re f ined model for the description of data. 
iKoposov et al.l ([2008H presented a search methodology 
for Milky Way satellite galaxies in SDSS data through 
the computation of efficiency maps. Search for stellar 
concentrations using these maps suggest a luminosity dis- 
tribution that steadily rises following a power law up to 



M„ 



5. From there on, a flat distribution could not be 



d iscarded (IKoposov et aTl l2008D. 

iTollerud et al.l (|2008( 1 use completeness limits for the 
SDSS-DR5 to implement a correction for luminosity bias. 
Although a first order correction would produce an in- 
crease in the faint end of the luminosity function, the 
authors bring forward that this result is not well enough 
constr ained given available data. iTrentham fc Tullvl 
( 2002) study the faint end of the galaxy luminosity func- 
tion in five different local environments, from the Virgo 
cluster to NGC 1023 group. The authors derive an aver- 
aged luminosity distribution in the range —18. < M r < 
— 10 (Cousins R magnitude) and infer a faint end slope 
a ~ —1.2. In the particular case of NGC 1023, a more 
detailed study confirmed later that the faint end is con - 
sistent with a sha ll ow sl ope ([Trentham fc Tullvl [20091. 
iTullv fc Trenthaml (pOOl studied the NGC 5353 group 
and attributed a faint end slope a — —1.15 to the fact 
that this group is at an intermediate evolutionary age. 

Membership of individual galaxies through spectro- 
scopic measurements is inefficient in terms of observ- 
ing time, given the large fraction of background objects 
that have to be rejected. This results in few systems 
with derived luminosity function complete down to faint 
magnitudes. For this reason, a background subtraction 
technique is an efficient method to study, on a statis- 
tical basis, properties of the population of companion 
objects. Using deep mock catalogues constructed from a 
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Figure 1. Properties of primary galaxies in the total sample. This 
sample comprises 51710 galaxies brighter than M r = —20.5 in the 
redshift ran ge 0.03 to 0.1, satisfying the isolation criteria stated in 
Section [2~H 



numerical simulation, iValotto et al.l ([2001( 1 analyze sys- 
tematic effects in the determination of the galaxy lu- 
minosity function in clusters. Their results indicate a 
strong tendency to derive a rising faint end when clus- 
ters are selected without redshift information. This is 
due to projection effects, since many of the clusters 
selected in 2D hav e no significant counterpart in 3D. 
iMuiioz et al.l (|2009l) use Mock catalogues constructed 
using GALFORM ([Baughl 1200(1 semianalytic model of 
galaxy formation to study the reliability of the statisti- 
cal background subtraction method to recover the un- 
derlying observer-frame luminosity function of high red- 
shift (z ~ 1) cluster galaxies in the K s band. These 
authors find that the optimal response of the method in 
recovering the underlying galaxy luminosity function oc- 
curs when background corrected counts of faint galaxies 
are complemented with photometric redshifts of bright 
galaxies. They also show that the increase in the number 
of galaxy clusters that contribute to the compu t ations 
dramatically reduce stochastic errors. IChristleinl IpOOOh 
study luminosity functions for galaxies in loose groups 
and suggests that the ratio of dwarf to giant galaxies is 
continuously increasing from low to high mass groups. 
This is based on data from the Las Campanas Redshift 
Survey, where environment is estimated using the line- 
of-sight velocity dispersion of the host groups. Back- 
ground subtraction has been applied to sin gle clusters 
(| Andreon et all 120051: iBarkhouse et al.ll2007f) and to en- 
sembles of clusters ( Gonzalez et al . 2006^ 

The goal of this work is to obtain statistical proper- 
ties of satellite galaxies in the magnitude range —18.5 < 
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M r < —14.5, using SDSS photometric data. Previous 
studies of galaxy satellites concern mostly Milky Way 
and nearby galaxies, so that this work complements with 
a study of projected radial number density profiles, lu- 
minosity, and colour distributions for a statistically large 
sample of primaries within z = 0.1. The definitions of 
host samples and the adopted satellite selection criteria 
are presented in Section[2j The details of the background 
subtraction procedure are given in Section [3] Then, in 
Section 2] we present the derived distributions of satellite 
properties. In order to validate the implemented method, 
we test it in Section [5] Finally, the main conclusions are 
provided in Section [5] 

2. PHOTOMETRIC AND SPECTROSCOPIC DATA 

We use the large database provided by the Sloan col- 
laboration (SDSS, iStoughton et al.l l2002t) . which pro- 
vides photometric information of objects down to faint 
magnitudes. This survey has been carried out using a 
dedicated 2.5m telescope (jGunn et~aL|[200l) . and com- 
prises digital photometric information of stars and galax- 
ies in 5 bands (jFukugita et al.lll996t iSmith et all 2002) 
reduced by an automated pipeline (|Lupton et al.l 12001 . 
The limiting magnitude is 22.2 in the r-band. A set 
of the brightest and more concentrated galaxies in the 
main gala xy sample has been s elected for spectroscopic 
follow up (Bl anton et al.ll2003al) . This leads to the spec- 
troscopic galaxy catalogue, which c ontains galaxies wit h 
Petrosian magnitude r < 17.77 (|Strauss et alj 120021 ). 
Both photometric and spectroscopic data are accessi- 
ble through a web interface to the public releases of 
the survey. The Data Release 7 comprises nearly 18 
Tb of catalogued data for objects identified as galax- 
ies by t he automated data red uction pipeline (Data Re- 
lease 7. IAbazaiian et al.ll2009l ). The main galaxy sample 
comprises 691055 galaxies with 95 per cent complete- 
ness down to a limiting magnitude rn m = 17.77 in the 
r-band. The surface brightness limits are imposed by 
the instrument capabilities and the automated reduction 
pipeline, so that these data sets allow to retrieve informa- 
tion about galaxies above the surface brightness limit of 
the cata logue, Mnn < 24.5 m ag arcsec~ 2 in the Petrosian 
r-band ([Strauss et al.l [2002) . All galaxies in this sam- 
ple have redshift measurements and serve in this work 
as primary targets for the study of fainter galaxies, ac- 
cessed from a deeper photometric sample. We obtain 
statistical properties of faint satellites, most of them not 
present in the spectroscopic survey, associated to pri- 
maries with measured redshifts in the range 0.03 to 0.1 
To this aim, we use the New York Value Ad ded Galaxy 
Catalog (NYU-VAGC, iBlanton et all l2005h to extract 
photometric information of neighboring galaxies of pri- 
maries. This catalogue is based on the sixth release of 
the Sloan Digital Sky Survey, and covers 9583 deg 2 of 
sky distributed into a northern cap and three southern 
stripes. We notice that the sky coverage of the main 
galaxy sample in DR7 is approximately included in the 
NYU-VAGC DR6 area, so that the cross-correlation be- 
tween these two catalogues is suitable for the purpose 
of th is work. We used software HEALPix (Gorsk i et al.1 
2005) to optimize the data processing of the galaxy cat- 
alogues. 

2.1. Host samples and galaxy selection criteria 



We have considered primaries brighter than 
M p = —20.5 (r-band luminosities) in the redshift 
range 0.03 to 0.1, applying an isolation criterion in 
order to avoid high density environments such as pairs 
or groups of galaxies. The density contrast around 
galaxies fainter than M p = —20.5 is low and comparable 
to Poisson uncertainty. Since the method, detailed in 
Section [3j is based on the presence of a strong signal 
to background, we left them out of consideration in 
defining the samples of hosts. These primaries are 
constrained to have no neighbors brighter than M p + 2 
within projected distance 700 kpc and relative radial 
velocity difference 700 km/s. These criteria are similar 
to those adopt ed in previous studies using spectro- 
scopic samples dSales fc Lam bas 2005]; iChen et all 120061 : 
lAgustsson fc Brainerdl2010h . and are intended to select 
and isolate halos where a dominating, assumed central 
galaxy of the satellite system is found. The total sample 
of primaries comprises 51710 objects brighter than 
M p = -20.5 in the redshift range 0.03 to 0.1. 

Taking into account galaxy luminosities and colours, 
we have considered different subsamples of primaries in 
order to explore possible dependencies of the satellite 
properties as a function of host properties. The descrip- 
tion of the subsamples considered is given in Table [TJ 
Subsample names are indexed on a three character basis: 
a number indicating the host luminosity selection ("0" 
for full sample, "1" for hosts with -21.5 < M r < -20.5, 
and "2" for hosts with M r < —21.5); and two letters 
indicating host and satellite colours. For simplicity, up- 
percase characters correspond to primary galaxy colours 
(A,R,B: All, Red, Blue), while lowercase indicate satel- 
lite colours (a, r, b standing for all, red and blue respec- 
tively). In Fig. [1] we show the distributions of r-band 
absolute magnitudes, g — r colours and redshifts of hosts 
in the total sample. The bimodal distribution of primary 
galaxy colours can be clearly appreciated, and we use a 
colour cut g — r=0.8 to separate subsamples according 
to host colour. 

3. BACKGROUND SUBTRACTION METHOD 

The background subtraction methodology are based on 
the simple idea of counting the number of objects in a 
region where a given signal is expected to lie, superposed 
to an uncorrelated noise, and subtracting a statistical es- 
timation of that noise. In this case, the signal is due to 
the presence of satellite galaxies in systems dominated 
by a central and luminous galaxy, and the noise is as- 
sociated to the background and foreground galaxies not 
dynamically linked to the primary galaxy. Then, it al- 
lows to statistically obtain properties of the faint galaxies 
associated to the primaries, without the need of redshift 
information for individual objects. This is accomplished 
provided that the working hypothesis of the central pri- 
mary is satisfied, and convenient ranges of observed pa- 
rameters are chosen, so that to minimize the contribu- 
tion of background counts. Although the method does 
not allow to quantify the contribution to the signal of 
each individual object, it is possible to obtain statistical 
estimates of probability density distributions describing 
galaxy properties. This is accomplished by restricting 
galaxies to a fixed bin of the variable, for instance the 
luminosity of the satellite or their distance to the central 
galaxy, and normalizing to the number of contributing 
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Table 1 

Definition of satellite samples considering colour range, maximum projected distance to the primary and primary properties. The number 
of primaries N p of each sample is also indicated. Satellite galaxies have a lower bound of g — r = —0.4 for samples a and b. Upper case 
letters A, B and R stand for All, Blue and Red, respectively. Similarly, lower case letters indicate the ranges of satellite colours. 



systems in that bin, after implementing the background 
subtraction. 

Within the hierarchical clustering paradigm, the mat- 
ter distribution can be roughly described by a set of halos 
populated with galaxies, according to certain reci pes that 
depend on galaxy type (jCoorav &; Shethl 120021 ) . This 
model can give insight to adopt an appropriate choice 
of the centers, which is key to increase the signal from 
the satellite population against the background noise 
and obtain reliable results using this method. Central 
galaxies play a different role than the rest of objects 
within the halo due to its particular accretion history, 
related to the merging of massive galaxies in each halo, 
and eventually the accretion of most of the available gas 
and e ven other minor galaxies (|Coorav fc Milosavlievid 

2005) . While galactic cannibalism has been proposed 
as the main mechanism for building up central galaxies 
(IQstriker fc Tremainell975llWhitel[l97alVale fc Ostriken 

2006) , it has also been suggested th at major merg- 
ers (ILin et"all 120041) and dry mergers (|Liu et al.1 120091 : 
iKhochfar fc Silkll2009f ) play an important role in the dif- 
ferent stages of their evolution. From a dynamic point of 
view, relative velocities of central galaxies with respect to 
the halos in which they reside are very small, compared 
to the large velocity dispersion of satellite galaxies. This 
leads to a clear observ ational distinction betw een central 
and satellite galaxies (jVale fc Qstrikerir2006| ). It is also 
known that a central gala xy is often the most lu minous 
galaxy within their halo (|Vale fc Ostrikerl 120061 ). Even 
when this " ce ntral galaxy paradig m" has been claimed to 
be inaccurate (jSkibba et al.l l2010). specially in high mass 
systems, the precise location of the center of the halo is 
not crucial for the implementation of the method, since 
it integrates galaxy counts in a region where overdensity 
signal is clearly present. Still, the location of the bright- 
est galaxy in this special type of galaxy systems, where 
a galaxy strongly dominates in luminosity and the total 
mass of the systems is quite low, is a very good approxi- 
mation to the position of the center of the galaxy system. 
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Figure 2. Colour distributions of faint galaxies with measured 
redshifts in SDSS, in three different absolute magnitude intervals. 
Curves show spline interpolation. 

Under these assumptions, the brightes t galaxies are ex- 
pected to reside in the centers of haloes pones fc Formanl 
ll984HLin et"aIll2004HSmith et al.ll2005l) . and can be used 
statistically to analyze galaxy overdensities associated to 
satellites. Accordingly, we limit the sample of centers 
to bright isolated galaxies, since they are likely central 
galaxies of small haloes. 

3.1. Constraints imposed by the data 

The main requirement for the success of a statisti- 
cal background subtraction method relies on a signif- 
icant mean overdensity around a given sample of pri- 
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Figure 3. Colour distributions of excess counts of companions, 
in three different magnitude intervals, of bright primaries (sample 
S2-A) resulting from a background subtraction calculation. 



maries. This fact leaded us to consider bright galaxies, 
M p < —20.5, since the neighborhood of fainter primaries 
do not show a significant density enhancement in SDSS 
data. The uncertainty of the background decontamina- 
tion procedure is dominated by Poissonian statistics of 
subtraction of large numbers. A significant increase in 
the signal to noise ratio of satellites can be achieved by 
eliminating those galaxies with a low probability of as- 
sociation to the primaries. To this end, we restrict the 
parameter space of galaxies in the photometric catalogue 
so that we only use ranges of those parameters where 
the hypothesis assumed in the background subtraction 
method are best satisfied, and reliable estimations of lu- 
minosity, colour and projected radial distributions can 
be obtained. We impose constraints on apparent magni- 
tude, colours and radial distance of photometric galaxies. 
Accordingly, we adopted M r < —14., —0.4 < g — r < 1.0, 
and a projected radial distance in the range 100 kpc to 
3 < R v i r >, which we discuss in detail in what follows. 

The redshift range of primaries determines, along with 
the limiting magnitude of the photometric sample, the 
maximum luminosity that can be studied. Given the 
limiting magnitude of 21.5 in the r-band, and a min- 
imum redshift for the primaries of 0.03, we chose the 
maximum luminosity M r — 5log(h) ~ —14 in the r-band. 
The limiting magnitude we use is lesser than the limit- 
ing magnitude from SDSS in order to ensure 100 per cent 
completeness (fAbazaiia n et al.ll2004[h 

We show in Fig. [2] the observed colour distributions 
of SDSS satellite galaxies derived from the spectroscopic 
sample in three different absolute magnitude bins be- 
tween — 18.5 and —14.5. It can be appreciated in this 
figure that a colour cut in g — r < 1 includes most of faint 
companion objects, and has the advantage of removing 
high redshift galaxies reddened by K-correction. There- 
fore, in all computations we exclude objects redder than 



this threshold to lower the noise due to the presence of 
high redshift galaxies. A clear indication that our colour 
cut g — r — 1 is a suitable threshold to remove high red- 
shift galaxies can be seen by inspection to Fig. [3] where 
we have applied background subtraction counts to sam- 
ple S2-A in the range 100 kpc< r p < 660 kpc. The ob- 
served lack of excess signal beyond g — r = 1 shows that 
our method is effective in detecting companion galax- 
ies in the colour range —0A<g — r<l. This is also 
a convenient cut according to determinations of satel- 
lite colours in sem ianalytic models of galaxy formation 
(iFont et all 120081) and galax ies in groups in SDSS-DR2 
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data (jWeinmann et al.ll2006D . We have applied the same 



colour cut g — r < 1 to galaxies in the mock catalogue, 
as described in Section [5] 

In Fig. |3]we show the projected density profile of galax- 
ies, in the selected color range —0.4 < g — r < 1.0, around 
primaries of sample S2-A-a. For comparison, we also 
show the corresponding of galaxies with g—r> 1. As can 
be appreciated, a significant excess of —0.4 < g — r < 1.0 
galaxies is observed beyond 100 kpc, while red galax- 
ies, g — r > 1.0 (mainly consisting in strongly redshifted 
background galaxies) show a null flat profile consistent 
with a uniform radial distribution (shaded region in Fig. 
21), which also gives support to our choice for the colour 
range of satellites. However, as can be appreciated in 
the inset of Fig. |4] we notice that the projected density 
profile of red galaxies is not uniform in the inmost region 
around the luminous primary galaxies. This shows that 
the hypothesis of a uniform background fails within 100 
kpc for this data set. Hence, in the computation of the 
statistical distributions, namely luminosity and colours, 
galaxies are restricted to be at a projected radial dis- 
tance of at least 100 kpc. This radial distance is con- 
sistent with previous findings that indicate the presence 
of an extended stellar halo associated to luminous galax- 
ies. This issue could lead to a failure of SDSS automated 
pipeline in the detection of low surface brightness galax- 
ies, either satellites or foreground /background galaxies. 
INierenberg et all (| 2 1 If) study the spatial distribution of 
faint satellites at intermediate redshifts (0.1 < z < 0.8) 
using high resolution HST images of early type galax- 
ies taken from GOODS fields. The authors model the 
light profile of host galaxies in order to study the pop- 
ulation of faint satellites. This model gives the spatial 
distribution of satellites near the primaries as a combi- 
nation of a satellite population with a power law radial 
profile superimposed to an isotropic and homogeneous 
background population. The method proposed by the 
authors requires the analysis of images of all hosts, which 
is beyond the scope of this work. This stellar halo com- 
ponent extends up to approximately 100 kpc, according 
to statistical determ i nations in SP SS (Zib etti et al.ll200"4T 
IBergvall et al.H2010HTal fc van Dokkumll201lD . 

An additional problem is related to fluctuations of lu- 
minosity density the extended images of luminous galax- 
ies, for example in the external parts of large, late-type 
galaxies at low redshift. Since SDSS images are analyzed 
to identify galaxies by automatic algorithms, luminosity 
concentrations (which are part of the primary) are of- 
ten confused and classified as many faint fake galaxies. 
This problem is hard to address, and is part of proposed 
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Figure 4. Mean projected density profile of galaxies with g — r< 
1.0 around primaries with M r < —21.5 (sample S2). Shaded region 
show the projected density profile of red galaxies around the same 
primaries with g — r> 1.0 (mainly background objects), within 1-cr 
uncertainties. 
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Figure 5. Density profile of satellites in S2-A sample selected 
by satellite luminosity. Galaxies are selected in the colour range 
—0.4 < g — r < 1.0. Fainter satellites show a steeper and more con- 
centrated density profile. Brighter satellites show a greater spatial 
extension, but smaller excess. 



SDSS "research challenges" Q The deblending process 
resulting from the current SDSS-DR7 pipeline appears 
to induce an artificial neighbor excess, which strongly 
difficults the tr eatment of data on reg ions close to the pri- 
mary (see e.g. iTollerud et al.1 (|2011[ )). specially in late- 
type galaxies where rich luminosity patterns are present 
on the disk. We also show in the inset of Fig. @]the mean 
value of the Rqo distribution of primaries i n sam ple S2, 
indicated with an arrow. Also, iBrainerdl [2005) argue 
that satellites show a strong preference for being aligned 
with the host major axis on scales below 100 kpc, while 
on larger scales (250 < r p /kpc < 500) the satellite distri- 
bution is consistent with an isotropic distribution. Due 
to this anisotropy, satellites close and near the major axis 
of the primary could be confused to luminosity enhance- 
ments of the disk of the primary when using photometric 
data. Both fake galaxies and depleted background are 
likely to coexist, and can not be directly disentangled. 
This might be also complicated by other effects such as 
dust obscuration by primary disk, alignment of satellites 
with disks and segregation of satellit e pro perties with 
radial distance (|Cnerj|2008t lAnn et al.H2008[ ). 

The chosen minimum distance of 100 kpc co rresponds 
roughly to 0.5 - 0.6R vlr . ITollerud et"aT1 (|2011l ) find that 
12 per cent of Milky Way-like galaxies host an LMC-like 
satellite within 75 kpc projected distance, while 42 per 
cent lie within 250 kpc. In the sample with the lower 
mean virial radius (SI) this means that, for a typical 
primary galaxy, and an assumed power law profile of the 
satellite distribution outside 20 kpc (7 = 2., see Table 
(2), only about 20 per cent of satellites within 100 kpc 
are not included in the analysis of sample SO and ~ 14 
per cent for sample S2. 

The projected density profile shows an excess of galaxy 
counts per unit area, with respect to the local back- 
ground, which decreases beyond R p « 500 kpc, but de- 
pends on the sample of prima ries and satellites consid- 
ered. iSales fe; Lam bas (2005) analize simulations and 
find that the turn around radius of satellite galaxies is of 
order 3 R V i r . These findings are also consistent with the 
approximate location of the turnaro und radius accordin g 
to the simple secondary infall model iBertschingerl ([1985). 
Then, in the computation of the statistical distributions, 
galaxies are also restricted to be at a projected radial dis- 
tance from the primary of up to 3 times the mean virial 
radius in each sample. Using the masses of dark mat- 
ter haloes hosting primary galaxies in a mock catalogue 
(described in Section [5]), we find a variation in the distri- 
butions of virial radius of primary galaxies corresponding 
to samples SO, SI and S2. While the median virial ra- 
dius < R V i r > for sample SO is 160 kpc, it changes to 
157 kpc for sample SI and to 220 kpc for sample S2. We 
adopt a maximum radius for each sample (SO, SI and S2) 
Rmax = 3 < R vir >, i.e., R max = 480 kpc for sample SO, 
Rmax — 470 kpc for sample SI, and R max — 660 kpc for 
sample S2. These values are also consistent w ith the re- 
sults of lChen et ail (|2006l ) and lLiu et all (f20ToT) , who find 
that the fraction of interlopers remain low within 500fcpc 
projected distance. In Section[5]we use a mock catalogue 
to test the background subtraction method, and confirm 
that the choice of this maximum radius is convenient to 
maintain a low level of contamination by foreground and 



1 http:/ /cas. sdss.org/dr7/en/proj/challenges/hii/default. asp 
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Figure 6. Density profile of excess galaxies around primaries in 
samples SO-A-a, Sl-A-a, and S2-A-a. For simplicity, only the power 
law fit is shown for sample SO-A-a, along with the uncertainty 
amplitude in long dashed lines. 

background galaxies. 

In Fig. [5] we show the resulting profiles for different 
ranges in satellite luminosities, where it can be seen that 
fainter satellites arc more strongly concentrated. This 
maximum projected radius 3 < R V i r > is suitable to 
study satellite properties in the range —18.5 < M r < 
— 14.5, although brightest satellites can b e det ected be- 
yond this value. Ag ustsson fc Brain crd (2006) analyze 
simulations in a LCDM framework to study the locations 
of satellite galaxies. The authors define their samples of 
satellites as all galaxies around luminous central galax- 
ies, that have a difference in projected distance lesser 
than 500kpc, and a difference in radial velocity below 
500-1000 km/s, according to the selected sample. Pri- 
mary gal axies have a med ian host halo virial radius of 
275 kpc. lAnn et all ((200 8) analyze satellite galaxies in 
SDSS-DR5 in order to study the dependence of radial dis- 
tribution and environment of galaxy morphology. These 
works show agreement on the extension of the satellite 
populations up to ~ 500 kpc. 

In order to account for large scale angular fluctuations 
in the distribution of faint galaxies we have used a lo- 
cal background in all samples. The density profile is 
approximately constant beyond 1.5 Mpc (see Fig. |6]), 
therefore we chose the density in the projected radius 
interval 2000 — 3000 kpc as the reference level of back- 
ground galaxy counts. Although a global background 
could be used instead, we argue that a mean number 
density of galaxies obtained locally is better suited to ac- 
count for possible irregularities in the number density of 
background galaxies. We use the redshift of each host to 
compute the absolute magnitudes corresponding to the 
excess galaxy counts in different apparent magnitudes, 
by assuming that these true companion galaxies are at 
the same redshift than the primary. We compute a com- 
posite luminosity distribution by averaging over a sample 



of primaries with a given set of properties. The number 
counts per magnitude bin on the resulting composite lu- 
minosity distribution estimates are normalized according 
to the appropriate limits of magnitudes in order to assure 
completeness: 



•V, ^ A 



(1) 



where iVj is the number of galaxies within the i-th mag- 
nitude bin for the total ensemble, and Nij is the number 
of galaxies in that magnitude bin for the j-th primary. 
The normalization for each magnitude bin is given by the 
fraction of fields contributing within the completeness 
limits, Cj, and by the fractional area, Fj, determined by 
the mask. 

3.2. Detailed masks 

Since we expect a relatively small number of satellites 
compared to the total background number counts, sev- 
eral primaries should be used in order to add up the 
signal over the background noise. Since the number of 
satellites grows approximately linearly with the number 
of hosts and Poisson errors scale as the square root of the 
number of galaxies, a large number of fields is required 
in order to override the satellite population number over 
the background counts. In Table [1] it can be seen that 
all samples comprise thousands of primaries, and in par- 
ticular, the smallest sample (S2-B) contains nearly 5000 
host galaxies. Fields of galaxies are extracted from the 
NYU-VAGC around the primary galaxies. The survey 
area from which the data is drawn has not, however, a 
simple geometry, and the detailed structure of holes and 
borders is an additional problem for the method used in 
this work. While the external borders are relatively sim- 
ple, the survey has small disconnected holes of a variety 
of shapes. The origin of these holes relies on the presence 
of several bright objects which blanket the background 
galaxies, in particular the faint ones. The most impor- 
tant are stars in our galaxy, but we can also mention the 
trails of solar system objects or artificial satellites, cosmic 
rays, etc. Moreover, in certain cases, there is a regular 
pattern of holes, sometimes associated to an incomplete 
structure of stripes. Therefore, a mask must be built for 
each of the fields, in order to correctly account for the ef- 
fective areas used in the method. As an example we show 
in Fig. 7(a) a typical field where a hole in the projected 
distribution of galaxies can be seen. The construction of 
the mask is based on the ansatz that, on a coarse resolu- 
tion, faint galaxies provide a dense, nearly uniform back- 
ground so that the holes can be directly associated to the 
absence of faint objects in a patch of the sky. In order 
to identify regions without faint galaxies we use a Monte 
Carlo method. The method consists in determining dis- 
tances from the faint galaxies in the photometric sam- 
ple to a large number of uniformly distributed random 
points within the area considered. Those random points 
that are separated from its nearest neighbor galaxy at 
least by a percolation radius ro, are used to character- 
ize the holes, i.e., holes are detected as concentrations 
of random points percolated by using this criteria. We 
have adopted a percolation radius r equal to a fraction 
/ of the mean angular galaxy separation < D > of each 
field. The value of / was obtained by diluting several 
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Figure 7. In left panel we show an example of a field centered on a bright galaxy affected by a mask, where each point represents a galaxy 
in the photometric sample. Using the Monte Carlo procedure to construct a detailed mask, as described in Section l3.2l an angular fraction 
is eliminated (dark points) in the region where a hole is located. In this typical example, we use a 98.6 per cent of the area in the field. 
On the right panel, we show the distribution of the fraction of areas of the primary fields F after application of the masks. Fields with 
F < 0.9 are excluded from the analysis. 



realizations of complete and dense fields of galaxies, and 
retrieving in each case the ideal limiting value for the 
percolation radius. We calibrated the value of this frac- 
tion, and found a mean value of / = 0.17, with a slight 
variation as a function of the projected galaxy density in 
each field (|Laresll2009[ ). Given the complicated shapes of 
identified holes, we simply eliminate angular sectors that 
contain the r andom points inside holes, as can be appre- 
In this Figure we show in grey points 
axies used in the background subtrac- 



ciated in Fig. 7(a 
the locations oTga 



tion procedure, and in dark grey points the locations of 
galaxies in the same angular sector that the mask hole 
that were not considered in such analysis. We computed 
the distribution of the fractions of angular areas finally 
used in the analysis presented in t he fol lowing sections. 
We show this distribution in Fig. 7(b) [ where we have 
indicated the location of the fraction ol used area of the 
field shown in Fig 



7(a; 



7(b) 



It can be seen from Fig 
the high level of completeness of the data, with 9Ci per 
cent of the fields having less than 5 per cent of their areas 
lost by the application of the masks. Furthermore, there 
is a negligible amount of fields with less than 90 per cent 
angular completeness, so that these fields are excluded 
from the analysis. Notice that this method is able to 
identify detailed small scale features in the angular dis- 
tribution of photometric galaxies. Given the resolution 
required to account for small hole s, a mask constructed 
by u s ing standard methods (e.g. lHamilton fc Tegmarkl 
[2001 IGorski et all 12003 ISwanson et alJl2008|) would be 
computationally expensive for reaching the same accu- 
racy than the used methodology. Besides, this procedure 
allows to adapt the percolation radius for each field, in- 



stead of using a fixed resolution on a pixelization scheme, 
and makes the computation of areas straightforward by 
using angular bins. 

4. RESULTS: STATISTICAL PROPERTIES OF SATELLITE 
GALAXIES 

4.1. Determination of projected radial distributions 

Although it would be desirable to consider satellites 
as close as possible to the primary galaxies, there are 
systematic detection biases that strongly limit this pos- 
sibility. We have also observed that close to primaries, 
the detection of faint objects can be strongly biased 
due to several facts such as obscuration, confusion, 
etc. Galaxies behind bright galaxy discs could also be 
covered by intrinsic absorption, producing significative 
changes in the observed magnitudes. Taking these is- 
sues into account, we have adopted a minimum dis- 
tance of 100 kpc to primaries for our analysis assur- 
ing reliable and systematic-free samples of faint objects. 
We considered galaxies in projected radial distance bins 
from the corresponding primaries to compute the av- 
eraged density profile for the different samples. The 
radial bins are chosen so that all the resulting rings 
have the same area. This partition allows to explore 
the inner region with more detail. In Fig. [6] we show 
the density profile of samples Sl-A-a and S2-A-a corre- 
sponding to all primaries with —21.5 < M r < —20.5 and 
—22.0 < M r < —21.5 respectively, without restrictions in 
the colours of primaries nor satellites. It can be appre- 
ciated the smoothly declining mean density profiles, al- 
though the extent of the overdense region is larger for 
the brightest hosts. We have also computed the density 
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Figure 8. Cumulative luminosity distributions of excess galaxies 
around primaries in the three defined samples (dubbed SO, SI and 
S2). Lines represent the best Schechtcr function fits to the obtained 
distributions (see text). 

profiles for red (0.4 < g — r < 1.0) and blue (g — r < 0.4) 
companion galaxies of bright primaries with M r < —21.5, 
(samples S2-A-r and S2-A-b). It can be seen by inspec- 
tion to Fig. [5] that the two profiles have a similar spatial 
extent and radial density profiles. 

We have fitted power law functions to the radial den- 
sity profiles p(r p ) — Ar~ n . The resulting values of the 
profile slope 7 are given in Table [5] where it can be seen 
that the companion galaxies of primaries show in gen- 
eral a concentrated distribution. lAgustsson fc Brainerdl 
(2010) investigate the locations of the satellites of rel- 
atively isolated host galaxies in the Sloan Digital Sky 
Survey and the Millennium Run simulation. The au- 
thors find that the distribution of the satellites within 
500kpc of red, high mass hosts with low star formation, 
differ from the distribution of satellites of blue, low mass 
hosts with low star formation. 

4.2. Mean number of satellites and luminosity 
distributions 

We compute the mean number of satellites beyond 
Rmin = 100 kpc and up to a projected radial distance 
equivalent to R max = 3 < R V i r > for each sample, no 
extremely red galaxies (g — r > 1) where included, since 
both restrictions where applied in the computations as 
previously explained (Scction l3.1[) . The mean number of 
satellites (shown in Table < N s > in the magnitude 
range M r < —14.5 is estimated for all primary samples, 
as: 



1 1 



(2) 



where N p is the number of primaries in the sample, N? 
is the number of galaxies inside the inner ring of the j-th 

field, and N 3 out is the number of galaxies inside the chosen 
outer ring. The area of the background, corrected by the 



Figure 9. Cumulative luminosity distributions of excess galaxies 
around bright primaries (M r < —21.5). Luminosity distributions 
for red and blue hosts (with a colour cut in g — T = 1.0) are also 
shown. 

mask, is A times the area enclosed between the chosen 
radius of the signal inner region. The isolation condi- 
tion ensures that bright galaxies chosen as centers are 
somewhat separated from other luminous galaxies, and 
is intended for selecting a given group of galaxies with 
similar properties, but does not affect the calculation of 
the luminosity function, given that it includes satellites 
at most as bright as M s = M p + 2 The faintest possible 
in sample SO has an absolute magnitude M p = —20.5, so 
that no satellites with M s < —18.5 are expected to be 
found in the spectroscopic sample. If galaxies brighter 
than this limit are considered, the sample of satellites 
would not be complete, so we limit the study of satellites 
to the range —18.5 < r < —14.5. We assign uncertainty 
estimates to the mean number of satellites using Pois- 
son errors resulting from the number counts of galaxies 
within Rmax from the primaries and the expected num- 
ber of background galaxies: 



1 

N~p\ 



-N j t 

ft out 



(3) 



We find that the number of satellites depends on pri- 
mary luminosity. As can be seen in Table [2l sample SI 
has on typically 1 satellite within the range of color and 
magnitude adopted. This number increases to nearly 6 
satellites on average for sample S2, equally distributed 
between red and blue satel l ites. T his result is consistent 
with that of iMaccio et al.l (|2010[) . who study numerical 
simulations of Milky Way sized haloes as predicted by 
Cold Dark Matter based models of galaxy formation and 
find that the number of satellite galaxies increase with 
halo mass. On the other hand, the number of companion 
objects of primaries on the same luminosity interval has 
not a clear dependence with satellite g — r colour index. 
In all samples of primaries, the numbers of blue and red 
satellites are comparable. We notice that this is valid for 
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Seiniplc 


< N s > 


T 


S0-A-a 


1.5 ± 0.4 


-2.0±0.2 


SO-A-r 


0.7 ± 0.4 


-2.1±0.5 


SO-A-b 


0.9 ±0.2 


-1.8±0.4 


Sl-A-a 


0.9 ±0.5 


-2.4±0.3 


Sl-A-r 


0.2 ± 0.5 


-2.0±2.0 


Sl-A-b 


0.7 ± 0.2 


-1.9±0.5 


S2-A-a 


6.1 ± 1.4 


-1.5±0.4 


S2-A-r 


3.4 ± 1.3 


-1.5±0.7 


S2-A-b 


2.7 ± 0.6 


-1.6±0.3 


S2-R-a 


7.8 ± 2.2 


-1.1±0.6 


S2-R-r 


4.2 ± 2.0 


-0.8±0.7 


S2-R-b 


3.7 ±0.9 


-1.8±0.2 


S2-B-a 


4.7 ± 1.8 


-1.9±0.5 


S2-B-r 


2.5 ± 1.6 


-2.1±0.7 


S2-B-b 


2.1 ± 1.0 


-1.5±0.6 



Table 2 

Mean number of satellites for each sample, and fit parameters for 
their background subtraction derived projected density profile. 
The mean number of satellites per host < N s > was calculated in 
the i — band. Errors are calculated as described in Section [3] In all 
cases, primary redshifts lie in the range 0.03 < z < 0.10 and the 
indicated number of satellites include objects in the magnitude 
range —18.5 < M < —14.5, with the g — r colour index limited to 
the range —0.4 < g — g < 1.0, and at a projected distance 
between lOOfcpc and Rmax (see Table [TJ. 

the ranges of projected radius and magnitudes used in 
this study (Tabled]). 

In the computation of the luminosity distribution, a 
convenient normalization is set for each magnitude bin, 
so that the mean number of satellites per magnitude in- 
terval per primary is obtained. This was achieved by di- 
viding the number counts of excess galaxies by the num- 
ber of contributing prima ries in each magn itude bin. We 
performed Schechter fits (|Schechterill976f) to the differ- 
ential histograms of magnitude distributions, and com- 
puted the faint end slopes, given in Table [2] For samples 
SO-A-a, Sl-A-a and S2-A-a, the luminosity distributions 
are consistent with a Schechter function with a faint- 
end slope of — 1.3 ± 0.2, adoptin g a universal M* value 
of -20.44 (jBlanton et al.ll2003bri . 

The results shown in Table [2] show a better sig- 
nal to noise ratio as more primaries comprise the 
subsamples. Bright primaries {M r < —21.5) host on 
average 7 satellites while low luminosity primaries 
(—21.5 < M r < —20.5) host a mean of only 1 satellite 
galaxy per primary. In Fig. [5] we show the r-band lumi- 
nosity distribution of satellites around bright primaries 
(M r < —21.5) and primaries with intermediate luminosi- 
ties (—21.5 < M r < —20.5, sample SI). We also show in 
this Figure Schechter function fits computed using the 
universal value of M*, whith a faint end slope param- 
eter a = — 1.3 ± 0.2, obtained using a maximum like- 
lihood method. These luminosity functions indicate a 
lack of a dominant population of faint satellites, which 
would be reflected in much larger negative values of 
the a parameter. These results contr ast with those 
obtained for groups/cluste rs of galaxies (|Popesso et al.l 
l2005t iGonzalez et al.l 120061) , where the faint component 
contribution M r > —18. to a double Schechter fitting 
gives slope values as steep as a ~ —2. 
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Figure 10. Cumulative luminosity distributions of blue (0.0 < 
g - t < 0.4), red (0.4 < g-r < 1.0) and all (0.0 < g - r < 1.0) 
excess galaxies around bright primaries in sample S2-A. 
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Figure 11. Colour-magnitude distribution of excess counts of 
— 18 < M r < — 15 companions of bright primaries (sample S2-A) 
resulting from a background subtraction calculation. 

We have also analyzed the dependence of the results 
on primary colour index. As mentioned in Section [2.1l wc 
have used the threshold g — r = 0.8 to divide the samples 
of primaries. Similarly, we find a larger population of 
satellites associated to red hosts, with a slightly stepper 
luminosity distribution at the faint end (Fig. [9|. We 
have studied the radial density profiles of red and blue 
satellites around the different samples of primaries, find- 
ing that the system of blue satellites is in all cases, more 
extended than that of the red ones by approximately 30 
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per cent. Luminosity distributions of red and blue satel- 
lites are shown for the brightest hosts in Fig. [TOl 

We also computed a colour-magnitude diagram for the 
satellites obtained by means of applying the background 
subtraction method simultaneously to these two vari- 
ables. The results are shown in Fig. [11] for sample S2-A- 
a, where it can be appreciated the smooth extension of 
the spectroscopic data onto fainter objects obtained by 
our statistical approach. 

5. TESTING THE METHOD WITH NUMERICAL 
SIMULATIONS 

The success of a background subtraction method re- 
lies on the signal-to-background strength from satellites 
around bright galaxies, providing the overdensity en- 
hancement obtained in the stacking procedure. However, 
the superposition of large scale structures projected onto 
the sky could affect the uniformity of the background. In 
order to estimate the ability of the method to correctly 
reproduce the actual distributions of satellite galaxies in 
SDSS-DR7 data, we have tested it on a mock catalogue 
derived from a numerical simulation using similar con- 
ditions than that applied to the observations. We con- 
structed the mock catalogue within a 7r/2 sterradians 
light-co ne, based on a sem ianalytic model of galaxy for- 
mation dCroton et al.ll2006l) at redshift zero in the Millen- 
nium simulation (jSpringel et al.ll2005f ). which is publicly 
available from the German Virtual Observatory 0. This 
solid angle corresponds to approximately half the area 
covered by the spectroscopic DR7 catalogue of galaxies 
and the z = snapshot was replicated 8 times along the 
axes to achieve a suitable depth. 

Since the output of the semianalytic model includes 
magnitudes in the ugriz photometric system, the mock 
spectroscopic catalogue is obtained directly by selecting 
galaxies brighter than the limiting apparent magnitude of 
the Sloan spectroscopic galaxy catalogue, r=17.77. From 
the mock spectroscopic catalogue we extract a sample of 
mock primaries using similar criteria than in Section [5] 
which will be used as centers in the following analysis. 
For each galaxy, a redshift is assigned by placing a fidu- 
cial observer at one corner, and determining the comov- 
ing distance to the observer and the peculiar velocity of 
the galaxy. The evolution corrected r-band magnitude 
is: 

M r = -2.5 log(i) + E - 5 log(», 

where E is given bv lBlanton et all ()2003bf ). 

For the adopted r-band limiting magnitude 21.5 of the 
photometric sample, the maximum redshift of the mock 
should be 1.2, corresponding approximately to the dis- 
tance at which an intrinsically luminous galaxy is ob- 
servable within the absolute magnitude range explored, 
— 18.5 < M r < —14.5. Although the lack of evolution in 
both galaxies and structure is a drawback of the mock 
catalogue, it serves as a strong test of the method given 
that in this case there is a larger clustering amplitude of 
high redshift structures (i.e. more structures along the 
line of sight) compared to a mock catalogue with consis- 
tent evolution in clustering. 

Since we have adopted a colour cut g — r < 1 in SDSS 
data to reduce background noise, we have performed an 

2 http://www.g-vo.org 



appropriate noise reduction in the mock catalogue using 
a Monte Carlo procedure in order to reproduce the con- 
ditions of the observations. We have considere d photo- 
metric redshifts obtained bv 10 'Mill et al.1 (|2010h (private 
communication) for the SDSS-DR7 galaxy sample to de- 
rive the fraction P(z) of galaxies with observed colour 
index g — r < 1 as a function of redshift. A suitable fit 
to this probability is given by P{z) = 1 - 2.381 (z - 0.08) 
which rejects about half the galaxies with g — r > 1 
at z ~ 0.3. We have adopted this statistical procedure 
instead of a direct filtering of galaxies by the observed 
colour in the photometric mock catalogue, given that this 
would be model dependent, also requiring reliable K cor- 
rection for the semianalytic galaxies. We also checked 
that our results are not strongly dependent on the pre- 
cise assumed P{z), so that little modifications in the fit 
have minimum impact in the obtained radial and lumi- 
nosity distributions. 

We have firstly tested if the galaxy density profile 
around primaries derived by the background subtraction 
method is able to reproduce the actual projected 3D pro- 
file. For this aim we have computed the projected radial 
distribution of galaxies around centers reproducing sam- 
ple S2-A-a in the mock catalogue using the real space 
positions in order to test the reproducibility of the re- 
sults through the background subtraction method. We 
show the results obtained for this sample (S2-A-a) in the 
mock catalogue since it presents the stronger signal. The 
results described in this section, however, does not de- 
pend on the chosen subsample. The criteria to select 
satellites in the mock 3D catalogue, was chosen so that 
a galaxy from the semi-analytic output is considered a 
satellite if it is within 3R V i r of the host dark matter halo 
of the primary galaxy, in three-dimensional space. This 
is consistent with the criteria adopted in Section 14.21 for 
the selection of the region where the signal is studied in 
the observational samples. In Fig. [T3] we show the good 
agreement of the projected 3D and background subtrac- 
tion derived profiles, indicating that the method is effec- 
tive in recovering the true projected profile of companion 
galaxies. 

The luminosity distribution of satellite galaxies were 
computed using the halo membership and galaxy lumi- 
nosities in the mock catalogue, and compared to the lu- 
minosity distribution of galaxies obtained with the back- 
ground subtraction in the region 100 < R p /kpc < 660 
corresponding to this sample. The cumulative luminosity 
distribution obtained through the background subtrac- 
tion procedure reproduces the true underlying distribu- 
tion remarkably well, as can be appreciated in Fig. 1131 
We also display in the inset of this Figure the differential 
distributions, which also present a general consistency. 
Although a small difference is present between the two 
samples, that does not affect the shape of the luminos- 
ity distribution and the determination of the Schechter 
parameter, since this is due to a small difference in the 
first magnitude bin. The method also succeeds in re- 
producing the discontinuity at M r ks —17.6, due to the 
absolute limiting magnitude in the parent semianalytic 
galaxy catalogue. 

These tests discard the possibility that background 
structures affect the shape of the obtained luminosity 
distributions, so that the faint end slopes are computed 
on firm statistical basis. We stress the fact, however, 
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that this procedure is reliable provided that adequate se- 
lection criteria have been imposed to the data. Since pri- 
mary galaxies in our samples span a redshift range much 
smaller than most of galaxies contributing to background 
counts, the procedure used to compute radial and lumi- 
nosity distributions performs in convenient conditions, as 
our tests in the mock catalogue have shown. This result 
is in conformity with previous tests of the method in less 
favorable conditions as in the determinat i on of cluster 
lumin osity functions {Valotto et al. 200l|; iMufioz et al.l 

USE!. 

6. DISCUSSION AND CONCLUSIONS 

We have carried out different statistical analyses to in- 
fer properties of faint satellite galaxies in the projected 
distance range 100 < r p /kpc < 3 < R V i r >, associated to 
bright primaries taken from SDSS with redshift z < 0.1 
To this end, we have implemented a background sub- 
traction method on faint galaxies with photometric in- 
formation, that are close in projection to galaxies with 
measured redshifts. The innermost region of the satel- 
lite systems is not accessible using the proposed back- 
ground subtraction method and data from the photo- 
metric galaxy catalogue. However, assuming a power law 
profile for satellites and considering objects outside 100 
kpc radial distance from the host, we can study more 
than 80 per cent of satellites. We have used a mock 
galaxy catalogue based on a semianalytic model of galaxy 
formation from the Millennium simulation to test the 
method with similar conditions than in the observational 
data. According to the results of the tests performed, the 
method is able to provide a good estimation of the true 
distributions of luminosities and projected radial galaxy 
density as a function of the distances to the host (Fig. [12] 
and Fig. fT3]) . In our mock catalogue we test how the pro- 
jection of background structures affect our measurement 
in a worst case scenario. We conclude that these struc- 
tures does not affect the shape of the luminosity distribu- 
tion, provided that central galaxies are sufficiently bright 
(M p < —20.5) and isolated, and a colour cut g — r< 1 is 
imposed to satellites. We also find it important to define 
an adequate maximum radius, using theoretical insight 
and mock catalogues to calibrate the values of R ma x for 
each sample. 

In all samples of primaries (defined in Table [lj we de- 
tect an excess of faint galaxy counts and we can deter- 
mine statistically the properties of companion objects 
associated to the central galaxy. We find that the ra- 
dial density profiles of satellites are consistent with power 
laws of the form p(r p ) = Ar^ 1 , with —2.4 < 7 < 1.1 and 
that the maximum extent and amplitude of the overden- 
sity depends on the primary luminosity and colour (see 
Fig. [6]), as well as on galaxy luminosity (Fig. [5]). The 
dependence of the number of satellites with M r < —14.5 
on host luminosity is strong: bright primaries with M r < 
—21.5 host on average approximately 6 satellites, which 
is reduce d to ~ 1 satellite for SI primaries (Table [2]). 
iLiu et ail (|2010( ) investigate the probabilities of finding 
a Milky Way like galaxy to host satellites with luminosi- 
ties similar to the Magellanic Clouds. The authors report 
that 8 per cent of these galaxies have at least 2 satellites 
similar to the LMC and SMC. In order to compare with 
these results we computed the mean number of satel- 
lites in the magnitude range —19.5 < M r < —14.5 for 
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Figure 12. Projected density profile of galaxies around primaries 
in the mock catalogue. Open circles correspond to results obtained 
by the background subtraction method, and filled circles corre- 
spond to the projection of the 3D density profile within the halos 
of each primary. We show in the inset figure, the region we ac- 
tually use in this work. Since the mock catalogue is not affected 
by spurious photometric galaxy detections, we could also show the 
good agreement between both profiles at projected radius below 
100 kpc. 

Milky Way like primaries (-21.4 < M r < —21.), within 
z = 0.1, using the previously described methodology. We 
find a (3 ± 2) per cent probability to obtain systems sim- 
ilar to MW-LMC- SMC. Th i s is co nsistent with the re- 
sults presented in ILiu et al.l (|2010l ). since we exclude in 
the analysis the central l OOfcpc region. 

Recently, iWang et al.l (pOlOh use a deep photometric 
sample around spectroscopically identified galaxies, and 
found that projected density profiles show a similar slope 
to the correlation function slope, independently of galaxy 
luminosity. Due to the different luminosit y and red- 
shift r anges considered between the work of IWang et al.l 
pOll and ours, a direct comparison is diffic ult to per- 
form. While this paper was being reviewed, IGuo et all 
(|2011h presented a complementary analysis of the lumi- 
nosities of satellites of SDSS primary galaxies, finding a 
general agreement with our results. 

The redshift range of the spectroscopic sample and 
the apparent magnitude limit of the photometric cat- 
alogue allows to obtain luminosity distributions in the 
range —18.5 < M r < —14.5. The derived luminos- 
ity distributions can be well described by Schechter 
function fits (Fig. [5J. Our findings indicate that 
faint end slopes of the satellite luminosity functions are 
slightly rising (a = —1.3 ±0.2). This is in agree- 
ment with the luminosity distributions of gala xies in 
the l ocal group, in the same magnitude range (|Mateol 
1998). This result is valid for all samples and indi- 
cates that the population of satellites of bright isolated 
primaries are consistent with the nearly flat faint end 
slope of the gl obal luminosity function as derived for 
the SDSS data (jBlanton et alj|2003bl iBaldrv et al1l2005l; 
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Figure 13. Cumulative luminosity distribution of galaxies around 
primaries in the sample S2 of the mock catalogue. Open circles cor- 
respond to the results from the background subtraction method. 
Filled circles are the actual luminosity distribution of satellites 
within three times the virial radius of the primary galaxies. In the 
inset we show the differential luminosity distribution on a lumi- 
nosity range beyond the completeness limit of the mock catalogue, 
where the background subtraction retrieve a distribution consistent 
with zero galaxies for M r > —17. in the Millennium semi-analytic 
catalogue. 



IMontero-Dorta fe; Pradal 120090 . These findings contrast 
with the results obtained by similar methods in samples 
of clusters and groups where a sig nificantly steep function 
is obtained (a ~ —2 to —1.5. e.g. Ide Propris et al.lll995t 
iPopesso et all 120051: [Gonzalez et al.ll2006[ ). These results 
are expected, given the evidence from semianalytic mod- 
els of galaxy formation that suggest that the total mass 
of a dark matter halo determines the normalization and 
shape of the luminosity function (|Macci6 eTalll2OT0h . 
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